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ABSTRACT 

A higher dimensional analogue of the KP hierarchy is presented. Fundamental con- 
stituents of the theory are pseudo-differential operators with Moyal algebraic coeffi- 



^ \ cients. The new hierarchy can be interpreted as large- iV limit of multi-component 
(gl(A^) symmetric) KP hierarchies. Actually, two different hierarchies are con- 
structed. The first hierarchy consists of commuting flows and may be thought 
of as a straightforward extension of the ordinary and multi-component KP hier- 
archies. The second one is a hierarchy of noncommuting flows, and related to 
Moyal algebraic deformations of selfdual gravity. Both hierarchies turn out to pos- 
sess quasi-classical limit, replacing Moyal algebraic structures by Poisson algebraic 
structures. The language of W-infinity algebras provides a unified point of view to 
these results. 



1. Introduction 



The technique of Large- A^" hmit is a magic that produces new continuous dimensions 
out of finite matrix indices. According to observations since the early eighties [1] , 
large- limit of sl(A^) looks like the algebra of infinitesimal area-preserving (or 
symplectic) diffeomorphisms on a two dimensional surface S: 

lim sl(iV) ~ sdifr(S). (1.1) 

iV— >oo 

If central elements (scalar matrices) are added, the limit becomes a Poisson algebra, 
i.e., the Lie algebra realized by a Poisson bracket on the surface: 

lim gl(iV) ~ Poisson(E). (1.2) 

iV— »oo 

These intriguing observations are based on the existence of the so called "sine 
generators" in these matrix Lie algebras. The theory of sine generators further 
tells us [1] that this is just a special, rather singular case; in general, the limit can 
also become a Moyal algebra: 

lim gl(iV) ~ Moyal(E). (1.3) 

N—^0O 

Moyal algebras are a kind of "quantum deformation" of Poisson algebras, which is 
a unique deformation subject to a set of natural requirements [2]. 

These observations, although mathematically slightly problematical [3], have 
been a very useful point of view for understanding higher dimensional integrable 
systems. For instance, large- A?" limit of the sl(A^) Toda field theory becomes a three 
dimensional nonlinear system (the dispersionless Toda equation) [4]. Similarly, 
large-A^ limit of two dimensional nonlinear sigma models with the pure WZW 
Lagrangian coincides with selfdual gravity (the selfdual vacuum Einstein equation) 
[5]. These two models are related to Poisson algebras, and well known to be 
integrable. Recently, a Moyal algebraic analogue of selfdual gravity is discovered 
from the same standpoint [6] and shown to be integrable [7] . 
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Inspired by this point of view, we present in this paper a higher dimensional 
analogue of the KP hierarchy. The celebrated KP hierarchy is known to possess 
a gl(iV) symmetric version, the so called 'W-component KP hierarchy" [8]. Our 
new hierarchy may be thought of as a kind of large- A?^ limit of the A'^-component 
KP hierarchy. Actually, rather than literally considering such limit, we directly 
construct the hierarchy replacing gl(A^) by the Moyal algebra Moyal(M) on a sym- 
plectic manifold M. For simphcity, we mostly deal with the planar case (i.e., the 
case where M is a two dimensional symplectic manifold), but all results can be 
extended straightforward to more general cases. 

The new hierarchy (which we call the "nonabelian KP hierarchy with Moyal 
algebraic coefficients") turns out to possess several novel characteristics. Of par- 
ticular interest is the fact that there are two distinct types of hierarchies within 
this framework. The first type of hierarchy consists of a commuting set of flows 
like many other integrable hierarchies. This hierarchy may be thought of as naive 
large-A^ limit of the A^-component KP hierarchy. The second type, meanwhile, is 
a hierarchy of noncommuting flows. This noncommuting hierarchy has a reduc- 
tion to a hierarchy of integrable flows including Moyal algebraic deformations of 
selfdual gravity. These Moyal algebraic deformations are somewhat distinct from 
Strachan's deformation. 

Another large-A^ limit, i.e., Poisson algebraic limit, can be realizes as quasi- 
classical {h 0) limit of these Moyal algebraic hierarchies. Quasi-classical limit of 
the commuting hierarchy gives a higher dimensional extension of the dispersionless 
KP hierarchy [9] . Similarly, quasi-classical limit of the noncommuting hierarchy is 
related to hierarchies of integrable flows in ordinary selfdual gravity [10]. 

This paper is organized as follows. Sections 2 and 3 are intended to be a 
review of key ideas and tools in the theory of integrable hierarchies and Moyal 
algebras. The contents of these sections, in particular Section 2, are mostly model- 
independent. Sections 4 and 5 are specialized to the Moyal algebraic hierarchies. 
Relation to Moyal algebraic deformations of selfdual gravity is discussed in Section 
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6. Quasi-classical limit is treated in Section 7. Section 8 is a summary of results 
from the point of view of W-infinity algebras. 



2. General structure of integrable hierarchies 



We here present a general framework for dealing with various integrable hierarchies 
on an equal footing. A prototype is the standard treatment of the ordinary and 
multi-component KP hierarchies [8] , which is now reformulated in an abstract and 
model-independent way. 

2.1. Lax and zero-curvature representations 

Integrable hierarchies usually possesses two distinct representations — the Lax 
representation and the zero-curvature (or Zakharov-Shabat) representation. Alge- 
braic structures of these representations are governed by a Lie algebra (of operators 
or of matrices) with a direct sum decomposition into two subalgebras: 



Let ( )± denote the projection onto each component. 

Lax and zero-curvature representations consist of differential equations ("Lax 
equations" and "zero-curvature equations" ) for t-dependent elements and Bi of 
the above Lie algebras: 



where a ranges over a finite index set, i over another (finite or infinite) index set. 
Let us call L^-'s "Lax operators" and BiS "Zakharov-Shabat operators." The Lax 
repretation can be written 



G-. 




La = Lait) eg, Bi = Biit) e g+, 




dLa 



[Bi, La-], 




dti 
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where ti denotes a time variable of the "i-th" flow. The zero- curvature equations 
(or "Zakharov-Shabat equations") can similarly be written 



The Lax and Zakharov-Shabat operators should be selected so that these equations 
be consistent. 

Example (KP hierarchy): For the KP hierarchy, Q is comprised of pseudo- 
differential operators (of finite order) in one variable x with scalar coefficients, 



n>0 

(E ^nd^) _ - (E <_i = E (2-5) 



n<-l 



The Lax and zero-curvature representations are formulated in terms of a single Lax 
operator L of the form 

oo 

L^d, + J29n+id-'' (2.6) 

n=l 

and and an infinite number of Zakharov-Shabat operators B^, n = 1,2,... given 

by 



Example (multi-component KP hierarchy): The formulation of the A^-component 
KP hierarchy is very similar, but based on pseudo-differential operators with NxN 
matrix-valued coefficients. Besides an analogue L of the Lax operator in the one- 
component KP hierarchy, we need N additional Lax operators Ua, ck — 1, . . . ,N 
to formulate the Lax and zero-curvature representations. Flows are numbered by 
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a double index i = {n.a), (n = 1,2,..., a = 1,...,N) and generated by the 
Zakharov-Shabat operators 

Bna = {L^Ua)>Q ■ (2.8) 



2.2. Dressing operator 

Integrable hierarchies of the above type allows another fundamental representa- 
tion — the dressing operator representation. The dressing operator is a i-dependent 
group element 

W = W{t) e exp G- (2.9) 

of the subalgebra Q- with which the Lax operators L^- are expressed in the "dress- 
ing" form 

= WCaW-'^, (2.10) 

where are "undressed" operators that are independent of solutions in consider- 
ation. With a suitable choice of the dressing operator, the Lax and zero-curvature 
equations can be converted into a set of differential equations of the form 

dW 

-^^BiW-WGi, (2.11) 

where Gi are a commuting set of operators that are also independent of solutions 
of the hierarchy. The Zakharov-Shabat operators turn out to be written 

Bi^{WGiW-'^)^, (2.12) 

so that one can eliminate Bi from the above equations. This results in the equations 

dW 



dti 



{WGiW-^)_W. (2.13) 
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The last equations give a well defined set of flows on the group exp Q- generated 

by g-. 

Example (KP hierarchy): In the case of the one- and multi-component KP 
hierarchy, 1^ is a pseudo-differential operator of the form 

oo 

p^ = i + 5]wna-". (2.14) 

n=l 

Remark: In the traditional theory of integrable hierarchies, Ccr's arc usually 
t-independent. We shall, however, see in later sections that this is too restric- 
tive. Furthermore, we shall extend the present setting to a case where Gj's are 
noncommutative. 

2.3. Factorization relation 

The above flows on exp can be translated into a kind of "factorization 
relation" in the group exp Q generated by Q (also called a "Riemann-Hilbert 
problem"). The following formulation is inspired by Mulase's work on the KP 
hierarchy [11]. [Remark: In general, such a group exp Q might not exist in a 
mathematically rigorous sense, but some realization as a "formal group" is always 
available as Mulase illustrates for the case of the KP hierarchy. Because of this, 
we use this somewhat loose notation and dare to call it a "group."] 

A fundamental factorization relation can be written 

W(t)E(t)W(0)-^ = W(t) e exp g+, (2.15) 

where 1^(0) is the "initial value" of = W{t) at i = 0, W{t) is a ^-dependent 
element of exp and E{t) is given by 

E{t) = exp (^Y^tiGi^ . (2.16) 



7 



(To be more precise, one has to enlarge ^+ into a "formal completion" with respect 
to time variables [11].) The above relation can be rewritten 



E{t)W{0)-'^ = W{t)-'^W{t), (2.17) 

which looks more like a "factorization" or a "Riemann-Hilbert problem." The direct 
sum decomposition at the level of the Lie algebra, (2.1), yields the direct product 
decomposition 

exp Q ~ exp Q- x exp Q+ (2.18) 

at the group level. The factorizability is thus ensured at least in a neighborhood 
of t = 0. Eventually, the flows on the space of dressing operators are converted 
into the action of E{t) on the coset space 

exp Q- ~ exp Q/ exp ^+ 

(2.19) 

W ^ [W]^Wexpg+ 

Let us briefly show how to prove the equivalence of (2.11) and (2.15). (For 
more details, see Mulase's paper [11]). 

(2.15) =^ (2.11): Differentiating both hand sides of (2.15) and recalling (2.16), 
one can derive the relation 

f^n^)^W(t)G,)Em(or = ?^. (2.20) 



dti J dti 



Using the factorization relation once again, one can eliminate the initial data W{0) 
to obtain 

Since the right hand side gives an element of Q+, the Q- component of the left 
hand side should be identically zero. This gives Eq. (2.11). 
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(2.11) =^ (2.15): Eq. (2.11) can be rewritten 

^^{W{t)E{t)) = Bi{t)W{t)E{t), (2.22) 

and repeated use of this formula yields a similar expression for higher derivatives, 

d d 



iW{t)E{t)) = Bi,...i,{t)Wit)E{t), (2.23) 



where Bi^...i^. lies in Qj^. Evaluated at t = 0, they give Taylor coefficients of 
W{t)E{t), 

W{t)E{t) =VF(0) + ^tiBi{{))W{Q) + • • • 

i 

l + ^iiSi(0) + ---j W^(O), (2.24) 

and the ffist factor 1 + XI tiBi{0) + ■ ■ ■ on the right hand side becomes an element 
of exp ^_|_. Writing this factor W{t), one finds that the last relation is nothing but 
the factorization relation. 

3. Pseudo-differential operators with Moyal algebraic coefficients 

Basic constituents of our new hierarchies, too, are pseudo-differential operators. 
Pseudo-differential operators in the A^-component KP hierarchy have NxN matrix- 
valued coefficients. We replace these coefficients by elements of a Moyal algebra 
Moyal (M). For simplicity, we here deal with the simplest case where M is a two 
dimensional plane. 
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3.1. MOYAL BRACKET AND STAR PRODUCT 



Let {y, z) be a pair of canonical coordinates on a two dimensional planar sym- 
plcctic manifold. The Moyal algebra consists of functions (or of formal power 
series) of (y, z) equipped with the Moyal bracket 



{a, h}n = ^ sinh 



2 \dydz' dzdy' ) 



a{y,z)b{y',z" 



(3.1) 



y'=y,z'=z 



where ^ is a nonzero parameter ("Planck constant"). Actually, we rather wish 
to consider the Moyal algebra of formal power series of {y,z) (and of fi). This is 
indeed possible, because the right hand side of the above formula is still meaningful 
in such a case. The "Planck constant" h then plays only a formal role. [Remark: 
The above definition of the Moyal bracket is somewhat unusual; usually, "sinh" is 
replaced by "sin" . This is just for simplifying notations, and one can go back to 
the ordinary situation by replacing h — > ih.] 

The Moyal bracket is a kind of "quantum deformation" of the Poisson bracket 



{a,b} 



da db da db 
dy dz dz dy ' 



(3.2) 



reducing to the latter in the classical limit. 



{a,b}n^{a,b} (h^O). 



(3.3) 



The canonical coordinates (y, z) are indeed canonical for both these Lie brackets: 



{y,z}fi^l, {y,z}^l. 



(3.4) 



Unlike the Poisson bracket, however, the Moyal bracket can be written as a (nor- 
malized) commutator 

{a,b}ji — ?i~^{a * b — b * a) (3.5) 
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of an associative product — the star product [1] [2] 

-h ( 92 92 



a*b = exp 



(3.6) 



2 \dydz' dzdy' 

The star product is nothing but the composition rule of Weyl-ordered operators. 

3.2. PSEUDO- DIFFERENTIAL OPERATORS WITH MOYAL ALGEBRAIC COEFFI- 
CIENTS 

Pseudo-differential operators with coefficients in the Moyal (or star product) 
algebra are, by definition, linear combinations 

m 

A= an&i (3.7) 

n=—oo 

of powers of with x-dependent coefficients taken from the Moyal algebra: 

an{h,x,y,z). (3.8) 

Multiplication of two pseudo-differential operators, which we write A* B and con- 
sider an extension of star product, is defined term-by-term as: 

(ad^) * (bd-,) = E Q « * (3-9) 

Actually, it is more natural to expand pseudo-differential operators in powers of 
Udx rather than dx- This indeed becomes crucial in the construction of integrable 
hierarchies, as we now show. 

We now specify a Lie algebra Q with direct sum decomposition into two sub- 
algebras for our new hierarchies. The Lie algebra ^, by definition, consists of 
pseudo-differential operators of the form 

m 

A^n-^ dnihdxr (3.10) 



11 



with Moyal algebraic coefficients that behave smoothly as — > 0: 

an = ai^\x,y,z) + Oih) (h^Q). (3.11) 

It is not hard to check that these pseudo-differential operators are indeed closed un- 
der star product commutator [^4, B] — A*B — B*A. The direct sum decomposition 
is basically the same as in the multi-component KP hierarchies: 

Q —Q>o © Q<-1, 
( )+ =( )>0 : projection onto {hda:f, {fidx)^ , 
( )- —{ )<-!'■ projection onto (^9a;)~^, (^5a;)~^, (3-12) 

The dressing operator should be a group element of Q<-i: 

W = exp^ h-^A, A e g<-i. (3.13) 

One can further expand W into negative powers of hdx, 

oo 

W^l + J2Mf^da^)~'', (3.14) 

n=l 

but the coefficients Wn will have singularities at h = that are rather hard to 
control. This is by no means a special property of the present situation, but 
common to this kind of formulation of /i-dependent integrable hierarchies. 



12 



4. Hierarchy of commuting flows 



We now proceed to the construction of the first hierarchy that consists of commut- 
ing flows. The construction is almost paraUel to the ordinary and multi-component 
KP hierarchies. 

4.1. Equations of dressing operator 

Let us first formulate the hierarchy in the language of the operator W with 
Moyal algebraic coefficients. The Lie algebra Q and its direct sum decomposition 
have been specifled in the end of the last section. As in the case of the A'"-component 
KP hierarchy, the index i of flows is a double index {n,a). Both components, 
however, range over all nonnegative integers as n = 0, 1, 2, . . ., a = 0, 1, 2, . . .. 
Generators Gi of flows are given by 

Gna = (^S.jV, (4.1) 

which obviously commute with each other. Differential equations of the dressing 
operator are given by 



dW 

'^'^na 

Bna^{W*Gna*W-')^^. (4.2) 

To see that these equations fall into the general framework of the previous sections, 
let us rewrite these equations as 

dW 1 1 

Otfia 

h-^Bna = {W * h-'Gna * W-^) . (4.3) 

Now hr^Gna is an element of Q. Since W is assumed to be in exp Q<^i, the 
conjugation W * h~^Gna * again belongs to Q. Therefore fT^Bna becomes 
an element of Q>q. Thus the consistency of the above equations is ensured. 
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These flows are highly nonUnear and nontrivial in general, but lowest flows 
contain trivial ones. For instance, the Zakharov-Shabat operators for (n, a) — 
(0,0), (1,0), (0,1) are given by 



Soo = l, Bio^hda^, Boi^y. (4.4) 



Therefore 



dW _ ^ dW_ _ dW dW_ _ dW_ 
dtoo ' 9^10 dx ' dtoi dz 

i.e., W is independent of too, and the flows of tio and toi are just translations in x 
and z. 

The factorization relation, too, can be formulated in the same form. The 
exponential operator E{t) is now given by 

E{t) = exp, h-' {J2 tnaihd.Yy'') , (4.6) 



where exp^ denotes the star exponential, 

oo ^ 

exp^ y4 = 1 + ^ — A (n-fold star product). (4.7) 

With this star exponential operator, the flows on the space exp G<-\ of dressing op- 
erators can be identifled with the action of E{t) on the coset space exp Q / exp Q>q. 

4.2. Lax and zero-curvature equations 
We now consider the operators 

L^W *hda:*W-'^, U *y*W-^ (4.8) 
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as analogues of the Lax operators L and Ua of the iV-component KP hierarchy. 
They can be written 



oo oo 

L^hd, + J2 9n+i im) -^ [/ = t/ + ^ K„ * L-^ (4.9) 

n=l n=l 

and the coefficients turn out to have smooth hmit as ^ — > 0, 

gn^9li\t,x,y,z) + 0{h), Un ^ u^n\t,x,y, z) + 0{h). (4.10) 

[Proof: The dressing operator W is assumed to be the star product exponential 
of an element of Q-i. Evaluating W * hdx * and W * y * by use of the 
general formula 

e-^Fe-^ = Y + [X,Y] + ^[X, [X, ¥]] + ■■■, (4.11) 

one can indeed prove the above fact. Q.E.D.] 
Zakharov-Shabat operators are given by 



S„« = (L-*C/")>o, (4.12) 

where we abbreviate — L * • ■ • * L (n-fold star product) and U°' — U * ■ ■ ■ *U 
(a-fold star product). 

As a consequence of (4.3), these operators turn out to satisfy the Lax equations 

QT QTJ 
h^ = [Bna,L], h^ = [Bna,U] (4.13) 

and the zero- curvature equations 

^^-^^ + [S^a,5n/3]=0 (4.14) 

with respect to star product commutator [A,B] — A*B — B*A. 
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4.3. Extended Lax equations 



We now introduce another pair of Lax operators. This is inspired by the 

treatment of quasi-classical limit of the KP hierarchy [12]. The extra Lax operators 
are given by 

M * E{t) * X * E{t)-^ * W-^, 

V *E{t)*z*E{t)-^ ^W-'^, (4.15) 

and satisfy the same Lax equations 

dM dV 

= [Bna, M], h-— = [Bna, V] (4.16) 

as well as the canonical commutation relations 

other commutators = 0. (4-17) 

Let us specify the structure of these extra Lax operators in more detail. By 
use of formula (4.11), one can calculate the "pre-dressing" by E{t): 

E{t) * X * E{t)-^ = ^tnaimr-^y'' + X, 
n,a 

E{t) * z * E{t)-' = Yl octnaihd^Ty''-^ + z. (4.18) 
n,a 

Dressing by W can easily be evaluated by the formulas 

W *x* = X + (element of Q<-i), 
W * z* = z + (element of G<~i), 

* (^i^a;) V * = * t^"- (4.19) 
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Thus, eventually, M and V turn out to be written 



oo 



n-1 



-n-1 



n,a 



n=l 



oo 



(4.20) 




and the coefficients hn and v^, like those of of L and C/, have smooth limit as 



Remark: Wc could have defined M = W * x * * z * W-\ but 

then extra terms emerge on the right hand side of the Lax equations. Exactly the 
same phenomenon takes place in the KP hierarchy [12]; the above definition of M 
and V mimics the formulation developed therein. In this respect, we should define 
L and U, too, via the same "pre-dressing" by E{t). This however gives the same 
result because hdx and y commute with E{t). 



By "noncommuting flows" we mean a set of flows with generators Gi that do not 
necessary commute with each other. In general, such noncommutativity causes 
serious troubles in formulating a Lax formalism. In the following case, however, the 
generators obey rather special commutation relations, so that we can write down 
Lax equations etc. in almost the same explicit way as in the case of commuting 
flows. 



K = h\i'{t,x,y,z) +0{h), 



Vn = vh {t,x,y,z) +0{h). 



(4.21) 



5. Hierarchy of noncommuting flows 
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5.1. Equation of dressing operator 

The second hierarchy consists of flows with three series of time variables, t = 
(ti, t2, . . .), p = {pi,P2, ■ ■ ■), and q = {qi, q2, ■ ■ ■)■ The flows are generated by the 
left action of the star product exponential 

OO 'OO oo 

E{t,p,q) = exp, h-'(J2tnm" - Y,PnZ{hd.X + Y.^ny{hd:,T) (5.1) 

n=l n=l n=l 

on the coset space exp Q/ exp Q>q. These flows can be converted into flows on 
the space exp G<-i of dressing operators by the factorization relation 

W{t,p, q) * q) * 1^(0, 0, 0)-i = W{t,p, q) G exp g>Q (5.2) 

connecting W — W{t,p, q) with the initial value at {t,p, q) — (0, 0, 0). The gener- 
ators of flows in E{t,p, q) are noncommutative, 

[n-\hd^y''y, h-\hd^rz] = n-\nd^y"+'', (5.3) 

but note that commutators are all central. Thus the algebraic structure is very 
similar to Heisenberg algebras. This is a key to the following calculations. 

Because of this mild noncommutativity, one can explicitly write down differ- 
ential equations satisfied by E{t,p,q) as: 

ind,:rE{t,p,q), 

^ 00 

[-Z(hd,r - I E ^rn{nd,r+-) * E{t,p, q), 

00 

{y{mr + ^Y.P'-(^^-r^") *E{t,p,q), (5.4) 



h 



n 



dE{t,p,q) 

dtn 
dE{t,p,q) 

dpn 

dE{t,p,q) 
dqn 
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[Proof: To derive these formulas, let us factorize E{t,p,q) in two different ways: 



oo oo 

Eit, p, q) = exp h-' (j2 tnihd.r - E Pn<f^d,r 

71=1 n=l 
^ oo oo 

m,n=l n=l 

00 oo 

= exp h-' (E tnimr + E invimr 

n=l n=l 

^ oo oo 

+ 2 E Pmqn{tld,r+^) *exph-'(^-J2PnZ{hd,r). (5.5) 

m,n=l n=l 

Here we have used the Campbell-Hausdorff formula 

e^e^ = exp(X + Y + hx,Y] + ■ ■ ■) (5.6) 

for star product. (Note that because of the above commutation relations of Heisenberg 
type, multiple commutators disappear.) The first factorized form of E{t,p,q) is 
suited for calculating p-derivatives, because p^s are included only in the first fac- 
tor, and all terms included therein commute with each other. The second formula 
of (5.4) can thus be proven. Similarly, the second factorized form of E{t,p, q) can 
be used to evaluate g-derivatives. For f-derivatives, nothing subtle is required. 
Q.E.D.] 

Having these differential equations for E{t,p,q), one can derive differential 
equations satisfied by the dressing operator as: 

dW 1 / i\ 

h^*w-^ ^- (w *(hda:T *w-^] 

Otn \ /<-l 

O T T 7" 1 OO 

* ly-i = - (-W * z(hdnT * w-^ - - y QmW * (hda:)'^+'' * w-^) 

Opn \ 2 ^-^^ /<-l 



dqn 



m=l 

oo 



m=l 
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[Proof: We simply repeat the reasoning in the end of Section 2. For instance, 
let us check the second equation. Differentiating the factorization relation and 
eliminating the initial value W{0, 0, 0), one obtains the basic relation 

(^ + W.^. E-A * ly-i = 1^ * W-\ (5.8) 

\dpn Opn ) Opn 

Since the right hand side should be an element of G>q-i 

(^left hand side of (5.8))^ ^ = 0. (5.9) 

By use of (5.4), one can easily check that this gives the second equation of (5.7). 
The other equations can be proven in the same manner. Q.E.D.] 

5.2. Lax and zero-curvature equations 

We now define Lax operators just hke the definition of M and V in the com- 
muting hierarchy: 

L * E{t,p, q) * hda: * E{t,p, q)'^ * W'^ , 
M * E{t,p, q)*x* E{t,p, q)~^ * W~^, 
U =W * E{t,p, q)*y* E{t,p, q)~^ * W~^, 

V =W * E{t,p, q)*z* E{t,p, q)-^ * W'^. (5.10) 
"Pre-dressing" by E{t,p, q) can be calculated by use of (4.11) as: 

E{t,p, q) * Wx * E{t,p, g)"^ = Wx, 

oo oo 

E{t,p,q)*x*E{t,p,q)-' = J^^tnikd^r-' - Y^np^zind^T-' 



n=l n=l 

oo 



+ Y^nqny{hdxr-'-^ J2 {m-n)pmqn{ridxr+^-' + x, 

n=l m,n=l 

00 

E{t,p, q)*y* E{t,p, q)'^ = ^Pn(^9^)" + V, 

n=l 

oo 

E{t,p, q)*z* E{t,p, q)-^ = Y, 9n(^5^)^ + z, (5.11) 



n=l 
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Dressing by W, then, yields the following expression of the Lax operators. 

oo 
n=l 

oo oo oo 

n=l n=l n=l 

, cx) oo 

+ 2 E {m-n)pmqnL"'+''-'+x + Y,hn+i*L-^-\ 

m,n=l n=l 
oo oo 

n=l n=l 

oo oo 

V^Yl + z + J2^n* L-^. (5.12) 

n=l n=l 

Just as in the commuting hierarchy, cf. (4.10) and (4.21), the coefficients on the 
right hand sides have smooth limit as /i — > 0. 

Let us now write down differential equations satisfied by these Lax operators. 
First note that (5.7) can now be written 

dW 



dtn \ )<-l 



dp^ 



n 



^ 00 



^dW 



m—l 

oo 



m=l 

Remarkably, M does not take place on the right hand side. This implies that one 
does not actually need M for defining Zakharov-Shabat operators. One can indeed 
derive, from the above equations, the Lax equations 



h— 

dtn 


= [Bn, 




dpn 






dqn 


^[Dn,Li 


^dU 

dtn 


= [Bn, 


u], 


^dU 

dpn 


— [Cm 


u], 


^dU 
dqn 


-[Dn,Ul 


^dV 

dtn 


= [Bn, 


V], 


^dV 

dpn 


= [Cn, 


V], 


^dV 
dqn 


= [Dn,V]. 



(5.14) 
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The Zakharov-Shabat operators Bn, Cn and Dn are given by 



m=l 

^ oo 



(5.15) 



m=l 



and satisfy the zero-curvature equations 



h 



h 



h 



h 



dBm 

dtn 

dBm 

dpn 

dBm 
dqn 
dCm 



n 



dpn 
dCm 



dqn 
dPm 
dqn 



- ^Trr^ + [Bm, Bn 

dtm 

- ^^TT^ + [Bm, Cn 

dtm 

^dDn p 

h [Bm, Dn 

dtm 

~ 1" [Cm, Cn 

dpm 

■ ^-K + [Cm, Dn 

OPm 

^dDn ^ r„ „ 

+ [Dm, Dn 

dqm 



=0, 

=0, 
=0, 
=0, 
=0. 



(5.16) 



Thus the Lax and zero-curvature equations are closed within L, U and V . (Of 
course M, too, satisfies the same Lax equations as the other three Lax operators.) 
Let us focus the following consideration on these three Lax operators. Note that 
they obey the commutation relations 



[L,c/] = [L,y] = o, [u,v\^n. 



(5.17) 
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6. Reduction to selfdual gravity 



Moyal algebraic deformations of selfdual gravity emerge from the hierarchy of 
{t, p, q) flows as a kind of "dimensional reduction" . Actually, connection to selfdual 

gravity takes two distinct forms, which correspond to two different local expres- 
sions of selfdual gravity — the first and second heavenly equations of Plebanski 
[13]. 

6.1. Constraints and reduced variables 

The reduction of the {t,p, q) flows is defined by the constraints 

^ = 0, n = l,2,... (6.1) 

on the coefficients of the dressing operator. This forces L and Bn to be trivial, 

L = nd:,, Bn = (^9^)". (6.2) 

In particular, 

dW dL dU dV , , 

^ = 0, — = — = — = 0. 6.3 

dtn Otn Otn Otn 

Under these constraints, hdx may be replaced by a "spectral parameter" A. The 
dressing operator then turns into a Laurent series with Moyal algebraic coefficients: 

oo 

W{X) = l + Y,Wn{h, t, p, q, y, z)X-^. (6.4) 
n=l 

This type of constraints are also used in the theory of multi-component KP hierar- 
chies [8] to derive an A'"-component version of the AKNS (Ablowitz-Kaup-Newell- 
Segur) or ZS (Zakharov-Shabat) hierarchy. The role of dressing operators is played 
therein by a Laurent series with matrix coefficients: 

oo 

^y(A) = l + ^u;^A-^ Wnit)eg\{N). (6.5) 

n=l 
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The Lax and Zakharov-Shabat operators, too, can be represented by Laurent 
series. As noted above, L is now trivial. Nontrivial dynamical contents are carried 
by U and V , which are now replaced by Laurent series of the form 

oo oo 

n=l n=l 

oo oo 

V{\) = ^ g,A- + ^ + J] ynX-"". (6.6) 



n=l n=l 



Algebraic relations among the Lax operators are carried over to these Laurent 
series. For instance, U (A) and V^(A) obeys the canonical commutation relation 

[U{\lV{\)] = h (6.7) 

and the dressing relations 

00 

U{\) ^W{\) * (j]pnA" + y) * iy(A)-i 

n=l 

oo 

-$^PnA" + iy(A)*|/*l^(A)-i, 

n=l 

oo 

V{X) =W{X) * (J2 QuX"" + ^) * W{X)-^ 

n=l 

oo 



n=l 



Laurent series corresponding to Zakharov-Shabat operators will be given in the 
next subsection. 
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6.2. Equations of flows in terms of Laurent series 



The factorization relation connecting W and its initial value is now translated 
into a factorization relation for 1^(A) = W{p, q, A) and its initial value aXp — q — Q: 

W{p, q, A) * Eip, q, A) * WiO, 0, Xy' = W{p, q, A), (6.9) 

where W{p, q, A) is a Laurent series to be obtained from W by replacing hdx A 
(which, too, turns out to be independent of x and t), and E{p, q, A) is the Laurent 
series 

oo oo 

E{p,q,X) =exp,h-^(^-J2PnZy' + Y.'l^y^'')- (6.10) 

n=l n=l 

(We put t — here and in the following.) The Lie algebra Q of pseudo-differential 
operators should thus be redefined to be a Lie algebra of Laurent series with Moyal 
algebraic coefficients: 

( )+ =( )>o : projection onto A^, A\ ■•■ , 

( )_ =( )<-! : projection onto A~^, A^^, • • • . (6-11) 



Differential equations for the Lax and dressing operators, too, can be translated 
into the language of Laurent series. First, W{X) satisfy the equations 



= " (t^(A)A'^)<-i * W{X). (6.12) 
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Similarly, U{X) and V{X) satisfy the Lax equations 



= i^niX), U{X)] = [- {V{X)X-)^, , U{X)l 



dpn 

dqn 

dp. 
dq. 

where Cn{X) and Dn{X) are given by 



= ^Dn{x\u{X)] = [{u{X)\-UM^)i 

= \Cn{X), V{X)] = [- (^(A)A")>o , ^(A)], 
= [Dn{Xl V{X)] = [(C/(A)A-)>o , V{X)l (6.13) 



Cn{X) = - {V{X)X-U + 9mA-+^ 

m=l 
^ oo 

DM = (C/(A)A")>o - - ^ p^A-+'^ (6.14) 



and obey the zero-curvature equations 



^aAn(A) _ ^9A^(A) ^ ^ ^ Q_ 

dqm 



6.3. Relation to second heavenly equation 

We now show that the above reduced hierarchy is related to a Moyal alge- 
braic deformation of self dual gravity. This deformation corresponds to Plebanski's 
"second heavenly equation" [13]: 

q2q q2q (dQ ^Q^ , , 



dydq dzdp \dy^ dz 
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Strachan's equation [6], on the other hand, is a deformation of Plebanski's "first 
heavenly equation" [13]: 



r on on r_ o'^n o'^n _ a'^n o'^n _ ^ 

I 9p ' dq J dpdp dqdq dpdq dqdp ' 



where {p, q, p, q) are another set of variables, and { , } stands for the Poisson 
bracket in {p,q). Comparing the above hierarchy with hierarchies constructed in 
Ref. 10, and identifying 

Pi=P, Qi=Q, (6.18) 

one will see that the above hierarchy gives a Moyal algebraic deformation of a 
hierarchy constructed therein for the second heavenly equation. Furthermore, Ple- 
banski's potential turns out to be linked with the next-to-leading term wi of 

w{xy. 

e = -hwi. (6.19) 

Let us show how to deduce these facts. To this end, we rewrite (6.12) as 

n^^B. * w{Xr' = - (^(A)A-)<_i , (6.20) 

and extract terms of both hand sides. This results in the relations 

dwi dwi 
h- — = Vn+i, n- — = -Un+i- (6.21) 
dpn Oqn 

[Actually, these relations can be extended to the case of n = by identifying 

Po = y, ?o = z, (6.22) 
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as one can directly check from the definition of U{\) and V^(A).] In terms of the 
function defined by (6.19), these relations can be rewritten 



Un = ^ , Vn = ■ (6.23) 

dqn-1 dpn-1 

In fact, exactly the same relations for selfdual gravity have been derived in Ref. 
10. Furthermore, repeating technical calculations presented therein, one can show 
that Eqs. (6.15) are actually equivalent to: 

dVm 9Vn . . „ 

+ {Vn,Vm}h = 0, 



dpn dpr, 

dUm dVn . . „ 

^ ^ ^{Vn:Um}n = 0, 

OPn Oqm 

dUm dUn 



dpn dpm 

By (6.23), these equations can be rewritten 



{un,Um}n^^, m,n = 0, 1,.... (6.24) 





J 


f de 


de 


dpndqm-l 


dpmdqn-1 1 


I dpn-1 


' dpm-l 


d^e 


d^e J 


f 90 


de 


dpndqm-l 


dqmdpn-1 1 


I dpn-1 


' dqm-1 




a^e J 




de 


dqndqm-1 


dqmdqn-l 1 


. dqn-1 


' dqm-1 , 



= 0, 
= 0, 

= 0. (6.25) 



The second equation with m — n — 1 gives a Moyal algebraic deformation of 
second heavenly equation (6.16). In the hmit of ^ — > 0, the Moyal bracket { , }fi 
turns into the Poisson bracket { , }, and these equations reproduce a hierarchy 
constructed in Ref. 10. 
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6.4. Relation to first heavenly equation 



Relation to the first heavenly equation is more involved. First heavenly equa- 
tion (6.17) is written in terms of the Poisson bracket { , } in (p, q), and Stra- 
chan's idea is to deform this bracket into a Moyal bracket as: 

I c?p ' dq J h dpdp dqdq dpdq dqdp 

These brackets differ from our brackets in {y,z). Nevertheless, very curiously 
and remarkably, it turns out that a similar Moyal algebraic analogue of the first 
heavenly equation emerges in a "gauge transformation" of the present hierarchy. 

Let us first specify this gauge transformation. By "gauge transformations" we 
mean symmetries of Eqs. (6.12-6.15) given by 



W{X) - 


W9{X) = 


= 9 


-^*W{X), 


U{X)- 


. U^{X) = 


9~ 


-^*U{X)*g, 


ViX)- 


V^iX) = 


9~ 


-'*V{X)*g, 


Cn{X) - 


- C^{X) = 


9~ 


* Cn{X) * g - g-'^ * 


Dn{X) - 


- D^X) = 


9' 


* Dn{X) * g - g~^ * 



dg_ 

dpn 

dg 

dqn 



(6.27) 



where g is a function of {y,z,p,q) and independent of A. With suitable choice of 
g (i.e., "a gauge-fixing"), A*^ terms in Cn{X) and Dn{X) can be eliminated. Such a 
gauge-fixing is realized by 

9 = ^0, (6.28) 

where wq is the leading term of the Laurent series 1^(A) that arises in the factor- 
ization relation: 

W(X)^wo + wiX^ + ---. (6.29) 
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The factorization relation, too, is gauge- covariant if W{X) is also transformed 



as 



W{\) W3{X) = g-' * W{X). (6.30) 

In the language of the factorization relation, the above special gauge transforma- 
tion simply means changing a normalization condition of solutions. By the above 
gauge transformation (note that g\p=q=Q — 1), indeed, the original normalization 
condition 

W^W|a=oo = 1 (6-31) 
is converted into a new condition of the form 

W^^(A)L=o = l- (6-32) 
This corresponds to taking a new direct sum decomposition of Q : 

Q —Q>i © Q<o, 
( )_|_ =( )>o : projection onto A^, A^, . . . , 
( )_ =( )<o : projection onto A°, A~^, . . . , (6.33) 

Repeating previous calculations in this new setting, one can derive the following 
equations, which are gauge transformations of (6.12). 



dpn \ ' ' J<0 



dW9{X) 



-(U3{X)X-)^^*W9{X), (6.34) 



dqn V y<0 

This kind of gauge transformations are also known for the KP hierarchy etc. [14] . 

Let us now return to the issue of selfdual gravity. We now augue that the 
function 

n = -hwf (6.35) 
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satisfies a Moyal algebraic analogue of the first heavenly equation, 

(mm\^ _d^^^_dH^d^ 0(h) -1 

\dp^ dq j dpdy dqdz dpdz dqdy 

where {pi,qi) are identified with {p,q) as in (6.18), and wf is the next-to- leading 
term of Laurent expansion of W^{X), 



W9(X) = i + wfX+---. (6.37) 

Let us show details. As a consequence of the factorization relation, W9{X) 
satisfies the same equation as W9{\): 



h ' = ( vs(x)x'') * W9(x), 

dpn V ' ' /<0 



Then imitating the derivation of (6.23), we can show that 



where Un and Vn are Laurent coefficients of U^^X) and ^^(A), 

oo oo 

C/5(A)=J]p„A" + «^ + J]<A-", 

n=l n=l 

oo oo 

V'{X) = J2 <ln>^'' + ^0 + E ^nA-"- (6.40) 

n=l n=l 

In particular, 

""0" dqi' '''dpi- ^^-^^^ 
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On the other hand, by the construction, Uq and Vq are hnked with the canonical 
conjugate pair y and z as 

Uq^9~^ *y* 9: v^^ g~'^*z* g, (6.42) 
hence canonical conjugate in themselves: 

K,^o'}n = l- (6.43) 

Substitution of (6.41) into this relation yields (6.36). Furthermore, as in the case 
of O, one can derive an infinite system of equations for Q. Those equations give a 
Moyal algebraic deformation of another hierarchy constructed in Ref. 10. 

We are thus in a somewhat puzzled situation. Our Moyal algebraic deformation 
(6.36) has exactly the same structure as Strachan's equation (6.26), but the Moyal 
brackets are living in apparently different planes, i.e., {y,z) and (p, g). In the 
standard treatment of selfdual gravity, the two coordinate systems (p, q, y, z) and 
(p, q) are distinct and never identical (unless the spacetime is flat). 

This problem will be resolved in a Toda version of our Moyal algebraic non- 
abelian KP hierarchies. The dressing operator approach of Ref. 7 , actually, em- 
ploys two dressing operators rather than a single one as we now considering. This 
situation resembles the Toda hierarchy [15], which, too, is based on two dressing 
operators. (Actually, although not clearly mentioned therein, fundamental ideas 
in Ref. 7 are rather borrowed from the theory of the Toda hierarchy.) A natural 
framework in which to embed Strachan's equation is thus a Moyal algebraic version 
of the Toda hierarchy. 
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7. Quasi-classical limit 



The Moyal algebraic hierarchies of both commuting and noncommuting types turn 
out to have quasi-classical {h 0) limit. Hierarchies in this limit, too, allow 
Lax- and zcro-curvaturc-typc representations, but commutators in the ordinary 
framework are now replaced by Poisson brackets. The situation is parallel to 
the KP hierarchy; its quasi-classical limit is the dispersionless KP hierarchy, and 
Poisson brackets are two dimensional. In the present setting, Poisson brackets are 
four dimensional. 

7.1. Prescription of quasi-classical limit 

Let us recall the derivation of the dispersionless KP hierarchy as quasi-classical 
limit of the ordinary KP hierarchy. The canonical conjugate Lax operators L and 
M of the KP hierarchy [12], in the presence of h, can be written 

oo 

00 oo 

M=^ni„L"-i + a; + ^/i„+iL-'^-\ (7.1) 

n=l n=l 

and satisfy the canonical commutation relation 



[L,M] = h (7.2) 

and the Lax equations 

^-[Bn,L], ^ = [S,,M], (7.3) 

where 
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The coefficients are assumed to behave smoothly as — > 0: 

9n = 9^n\t,x) + 0{n), K = hi^\t,x) + 0{h). (7.5) 

In the hmit of /i — > 0, the canonical conjugate pair {hdx,x) of operators are 
replaced by a canonical conjugate pair [k, x) of coordinates on a two dimensional 
symplectic manifold. The role of the Lax operators L and M are then played by 
the Laurent series 

oo 
n=l 

oo oo 

M=J2 ntnC""-^ + ^ + E (7-6) 

n=l n=l 

which satisfies the canonical Poisson relation 

{C,M}-1 (7.7) 

and the Lax-like equations 



with respect to the Poisson bracket 

fAB\-——-— — 

dk dx dx dk 

Here BnS are polynomials of the form 

and the projection operators ( )>o and ( )<-i are defined by 

( )>o : projection onto A;^, k^, - ■ ■ , 
( )<_i : projection onto k~^, k^^, ■ ■ ■ . (7-11) 
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These polynomials Bn then satisfy the zero-curvature-like equations 

dBm dBn 



dtn dtr 



+ {Bm,Bn} = 0. (7.12) 



Thus one can reproduce the Lax formalism of the dispersionless KP hierarchy. 

In the case of our higher dimensional hierarchies, we have anther canonical 

conjugate pair {y,z) that commute with {hdx,x), hence altogether four operators 
{hdx,x,y, z). These operators obey the canonical commutation relations 

[hdx,x]^h, [y,z]^h, 

other commutators = 0. (7-13) 

In the limit of ^ — > 0, these operators will be replaced by canonical coordinates 
{k, X, y, z) on a four dimensional symplectic manifold, and by the standard rule 

n-\ , ]^{{ , }} (7.14) 

of quantum mechanics, the commutator of operators will be reduced to the four 
dimensional Poisson bracket 

OAdB OA OB OA OB 
dk dx dx dk dy dz dz dy 

of functions of {k,x,y,z). The role of Lax and Zakharov-Shabat operators will 
thus be played by functions of {k,x,y,z). Those classical analogues of Lax and 
Zakharov-Shabat operators will satisfy Poisson algebraic analogues of the Lax and 
zero-curvature equations with respect to {{ , }}. Let us show details below. 
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7.2. Hierarchy of commuting flows 



For the hierarchy of tna flows, classical counterparts of L, M, U and V are 
given by the Laurent series 

00 

n=l 

oo 

,(0) 



n,a n=l 

oo 



n=l 

oo 

n,Q; n=l 

where the coefficients g'i^'', u^^ and f^'^'' arc the leading terms of h expansion 
in (4.10) and (4.21). These Laurent series satisfy the canonical Poisson relations 

{{A-M}} = {{W,V}} = 1, 

other Poisson brackets = 0, C^-l^) 
and the Poisson algebraic Lax equations 



f)f II- nut, ~jj, o, 

dU dV 



{{BnaM}}. ^ = {{Sna,V}}. (7.18) 

Here the classical counterparts Bna of the Zakharov-Shabat operators B^a are now 
given by 

B„« = (£-W")>o, (7.19) 

where 



)>0 : projection onto k ,k , , 



,-1 1.-2 



( )<_i : projection onto k ,k , • • • , (7.20) 
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and satisfy the Poisson algebraic zero- curvature equations 



r 



+ {{Bma,Bn0}}^O. (7.21) 



Note that if we now impose the constraints 

^ijf) = 0, v^^^ = for n = 1, 2, . . . , 
U = y, V = ^aW/:V"' + ^) (7.22) 

only the tno flows remain nontrivial. The reduced hierarchy is substantially the 
same as the dispersionless KP hierarchy. 

7.3. Hierarchy of noncommuting flows 

For the hierarchy of {t,p,q) flows, classical counterparts of L, U and V are 
given by 

oo 
n=l 

oo oo 

.(0), 



n=l n=l 

oo 00 

V=J2 + Z + ^n^^-"". (7.23) 



n=l n=l 

They satisfy the Poisson bracket relations 

{{C,U}} = {{C,V}} = 0, {{U,V}} = 1, (7.24) 
and the Poisson algebraic Lax equations 

rflJ BIJ rSlA 

%^^{{BnM}}. n^^{{CnM}}. h^^{{Vr,M]], 
dtn Opn Oqn 

n^ = {{Bn,V}}, f^^ = {{Cn,V}}, h^ = {{Vn,V}}. (7.25) 
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where Bn, Cn and V are given by 



BnHn>o^ 

oo 

m=l 

^ oo 

Vn^(uC---J2PmC'^+'')^^, (7.26) 



m=l 



and satisfy the Poisson algebraic zero-curvature equations 

+ {{Bm:Bn}}=0, 
+ {{Bm,Cn}}^0, 
+ {{Bm,Vn}} ^0, 
■ + {{Cm,Cn}}=0, 
+ {{Cm,'Dn}} =0, 



dBm 


dBn 


dtn 


dtm 


dBm 


dCn 


dpn 


dtm 


dBm 


dVn 


dqn 


dtm 


dCm 


dCn 


dpn 


dpm 


dCm 


dVn 


dqn 


dpm 


dVm 


dVn 



Selfdual gravity can now be reproduced under the constraints 

gl^l^^l forn = l,2,... = (7.28) 

All tn flows then become trivial, and upon identifying \ = k, the reduced hierarchy 
of remaining {p, q) flows reproduces the same hierarchy as constructed in Ref. 10 for 
the second heavenly equation. Note that in the previous section, wc flrst imposed 
constraints, then took the ^ ^ limit. In this section, we first took the ^ ^ 
limit, then put the constraints. We have thus two routes to the same second 
heavenly equation. 
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The status of the first heavenly equation is rather obscure from this point 
view. As mentioned in the last section, a better approach will be to start from a 
Toda version of the hierarchy of {t,p, q) flows, and to consider quasi- classical limit. 
Similarly, a Toda version of the hierarchy of tna flows will yield, in the same quasi- 
classical limit, a higher dimensional extension of the dispersionless Toda hierarchy 
[16]. 

7.4. MOYAL vs. POISSON ALGEBRAIC HIERARCHIES 

Remarkably, unhke the Moyal algebraic hierarchies, these Poisson algebraic 
hierarchies are systems of "algebraic differential equations" — they consist of an 
infinite number of equations, but each equation contains only a finite number of 
terms. This is not the case for the Moyal algebraic hierarchy, because the Moyal 
bracket itself is an infinite power series of h and contains arbitrarily high order 
derivatives of the unknown functions. First of all, Moyal algebraic deformations of 
the two heavenly equations themselves are non- algebraic in this sense. 

A drawback is that the dressing operator method does not work directly within 
the Poisson algebraic hierarchies. The dressing operator method can be used only 

via the Moyal algebraic hierarchies. Twistor theory, which has been successful for 
selfdual gravity and the dispersionless KP and Toda hierarchies [10] [12] [16], will 
provide an alternative approach. 

8. Conclusion — overview from W-infinity algebras 

Inspired by the philosophy of large- limit, we have constructed higher dimen- 
sional analogues of the KP hierarchy. Conceptually, the new hierarchies may be 
interpreted as large- A^ limit of the A^-component KP hierarchy. We, however, have 
avoided to deal with finite- A" models and given a direct construction of such hierar- 
chies. Fundamental constituents of the hierarchies are pseudo-differential operators 
with Moyal algebraic coefficients. 
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Actually, we have found two different types hierarchies consisting, respectively, 
of commuting flows and noncommuting flows. The hierarchy of commuting flows 
resembles the ordinary and multi-component KP hierarchies, whereas the hierarchy 
of noncommuting flows includes Moyal algebraic deformations of selfdual gravity 
as a reduction. 

Furthermore, both the commuting and noncommuting hierarchies have turned 
out to possess quasi-classical limit, whose Lax formalism is based on Poisson brack- 
ets rather than Moyal brackets. In the quasi-classical limit, the commuting hierar- 
chy gives a higher dimensional extension of the dispersionless KP hierarchy. The 
noncommuting hierarchy, meanwhile, includes ordinary selfdual gravity as a reduc- 
tion. 

It is instructive to review our results from the point of view of W-infinity 
algebras [17]. The following is a list of integrable systems that fall into our current 
scope. 

(0) ordinary KP hierarchy (KP), 

(1) A'"-component KP hierarchy (iV-KP), 

(2) nonabelian KP hierarchies with Moyal algebraic coefficients (MAKP), 

(3) Moyal algebraic deformations of selfdual gravity (MASDG) , 

(4) AT-component AKNS-ZS hierarchy (A^-AKNSZS). 

These integrable systems are linked with each other as the following diagram indi- 
cates. 

KP AT-KP ""-^ MAKP 



A^-AKNSZS MASDG 

[Remark: More precisely, large- A" hmit of A"-AKNSZS gives a reduction of the 
commuting MAKP hierarchy rather than of the noncommuting MAKP hierarchy. 
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hence the last hne of the above diagram is somewhat inaccurate.] The following list 
shows infinite dimensional Lie algebras associated with these integrable systems. 



(0) Woo = {x\nd:,y) (ordinary W^o algebra [17]), 

(1) — (Eai3X^{h,dxy) {Woo algebra with inner symmetries [18]), 

(2) W^ = {y'^z'^x'ihda^y) (large-A^ limit of W^ [18]), 

(3) >CMoyal(E) = (y^z^A^) (loop algebra of Moyal algebra), 

(4) >Cgl(A^) = (E^'f^X^) (loop algebra of gl(A^)), 

where < • • • > shows generators; E"'^ denote the standard generators of gl(A^), 

., = 6iaSjp. (8.1) 

The first four of these algebras are of W-infinity type, whereas the last one is of 
Kac- Moody type. The above diagram of integrable systems simultaneously show 
the interrelation of these Lie algebras: 

iV-comp. N-^ 



hdx 



jCg\{N) £Moyal(E) 

In the hmit of ^ ^ 0, these W-infinity algebras are contracted to classical counter- 
parts. In particular, £Moyal(S) and W^ turn into Poisson algebras £Poisson(S) 
and ~ Poisson(M), dimM = 4. This is exactly what we have seen in the 
transition from Moyal algebraic hierarchies to Poisson algebraic hierarchies. 

We conclude this paper with a list of issues to be pursued along the same line 
of approach. These issues will be discussed elsewhere. 

(a) Toda version of MAKP. We noted in a previous section that a Toda ver- 
sion of MAKP will give a natural framework for dealing with the first heavenly 
equation of selfdual gravity. This should be a rather straightforward extension, 
simply replacing pseudo-differential operators by difference operators (with Moyal 
algebraic coefficients). 
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(b) More detailed study of Poisson algebraic hierarchies. The Poisson alge- 
braic hierarchies will provide new material for the twistor theoretical approach to 
nonlinear integrable systems. Furthermore, as the dispersionless KP hierarchy is 
applied to Landau-Ginzburg models of topological strings [19], the higher dimen- 
sional hierarchies might be related to similar physical models. 

(c) Multi- component version of MAKP. This extension is also straightforward, 
just redefining the pseudo-differential operators W, L, etc. to have N x N matrix 
valued coefficients with Moyal algebraic matrix elements. According to a prelimi- 
nary analysis (we omit details), a hierarchy of Bogomolny-type [20] can be derived 
from this A^-component version (A^-MAKP) as a reduction. Thus the previous 
diagram can further be enlarged as: 

MAKP "^"-^P- AT -MAKP 

reduction 

Bogomolny 

This reduction is however somewhat strange — constraints are imposed on both 
the dressing operator and the time variables. This interpretation of Bogomolny hi- 
erarchies resembles the "D-module approach" of Sato et al. [21]. Furthermore, the 
fact that the time variables, too, have to be constrained is somewhat reminiscent 
of the "small phase space" of Landau-Ginzburg models [19]. 

(d) Hierarchies with compactified Moyal algebras. In the case of a torus, the 
Moyal algebra (or "quantum torus algebra") has an invariant trace [22] This will 
allow us to construct a higher dimensional tau function, because the usual trace 
of gl(A^) plays an imphcit but substantial role in the theory of tau functions of 
multi-component KP hierarchies. 

The author is very grateful to Jens Hoppe, Stanislav Pakuliak, Ian Strachan 
and Takashi Takebe for a number of useful comments. This work is partly sup- 
ported by the Grant-in-Aid for Scientific Research, the Ministry of Education, 
Science and Culture, Japan. 
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